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PURPOSE 


The  purpose  of  the  project  is  to  increase  understanding  of  the 
roles  of  cloud  conversion,  accretion,  evaporation,  and  entrainment 
processes  in  shaping  the  distributions  of  water  vapor,  cloud,  and 
precipitation  associated  with  tropical  circulations. 

ABSTRACT 

A  set  of  equations  describing  kinematical  relationships  among 
cloud,  precipitation,  and  the  wind  field  are  derived  in  detail. 
Steady-state  profiles  of  cloud  in  a  non-precipitating  tropical  at¬ 
mosphere  are  presented  for  10°C  steps  of  the  condensation  levelj 
the  unresponsiveness  of  cloud  content  to  vertical  air  displacements 
at  high  altitudes  and  low  temperatures  is  an  explanation  for  the 
relative  persistence  of  cirrus  clouds.  Exact  solutions  are  pre¬ 
sented  for  a  cloud  and  precipitation  model  that  is  simplified  to 
include  only  the  problem  of  precipitation  initiation.  Ibis  model 
shows  that  where  two  processes  for  precipitation  formation  exist, 
the  onset  of  precipitation  is  controlled  only  within  certain  limits 
by  alteration  of  one  of  the  processes.  Microphysical  parameters 
influence  transient  features  of  the  development  of  precipitation 
and  cloud.  In  many  cases,  however,  the  steady-state  surface  pre¬ 
cipitation  rate  is  only  slightly  responsive  to  order-of-magnitude 
changes  of  the  microphysical  parameters. 
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PUBLICATIONS .  LECTURE.  RETORTS.  AND  CONFERENCES 


lhe  paper  "Elementary  theory  of  association*  between  atmospheric 
motions  and  distributions  of  water  content"  by  E.  Kessler  was  pub¬ 
lished  in  the  January,  1963  issue  of  the  Monthly  Weather  Review. 

The  paper  "Relationships  between  tropical  precipitation  and  kine¬ 
matic  cloud  models"  by  E.  Kessler  was  presented  at  the  43rd  Annual 
Meeting  of  the  American  Meteorological  Society  in  New  Tork  City  on 
21  January  1963. 

Dr.  Kessler  and  Mr.  Feteris  met  with  Dr.  H.  Weickmann  at  the  New 
Tork  meetings  and  discussed  the  work  of  the  contract  and  future  plans. 

The  program  committee  of  the  10th  Weather  Radar  Conference  has 
accepted  our  paper  for  presentation  on  April  23,  1963.  The  abstract 
follows: 


Role  of  Microphysical  Processes  in  Shaping 
Vertical  Profiles  of  Precipitation  and  Cloud 

by 

Edwin  Kessler,  Pieter  J.  Feteris,  Edward  A.  Newburg 

Continuity  equations  are  presented  which  relate  cloud  and  preci¬ 
pitation  development  to  the  wind  field  and  processes  of  condensation, 
evaporation,  conversion  of  cloud  to  precipitation,  and  collection  of 
cloud  by  precipitation.  The  shape  in  steady  updrafts  of  transient 
features  of  cloud  and  precipitation  profiles  is  related  to  cloud  corn- 
version  and  collection  rates,  which  control  the  onset  of  precipita¬ 
tion.  However,  if  one  of  these  processes  is  relatively  rapid,  the 
rate  of  development  of  precipitation  is  quite  insensitive  to  the 
magnitude  of  the  other.  After  precipitation  starts,  hydrometeor  pro¬ 
files  approach  a  steady  state  in  all  cases  where  the  maximum  steady 
updraft  is  less  than  typical  precipitation  fall  speeds.  The  equations 
constitute  a  model  precipitation  system  whose  reaction  to  changes  of 
the  microphysical  processes  is  usually  in  a  direction  to  minimize  the 
effects  of  the  changes.  The  steady-state  vertical  profiles  are  in¬ 
dependent  of  initial  conditions,  weakly  dependent  on  the  microphysical 
parameters  and  strongly  dependent  on  the  updrafts. 
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1 .0  FACTUAL  DATA 


1.1  Comprehensive  Privation  of  the  Basic  Equation* 

A  detailed  derivation  of  the  continuity  equations  which  are  the 
foundation  of  this  study  has  not  been  given  previously.  This  deriva¬ 
tion  is  of  particular  importance  to  two-dimensional  problems  now  being 
designed  and  programed  for  solution  by  digital  computer. 

1.1.1  A  continuity  equation  for  precipitation 

Assume  that  the  fall  speed  of  the  precipitation  content  M 
relative  to  the  air  at  any  one  height  and  time  can  be  represented  by 
the  parameter  Vj  if  the  precipitation  content  at  a  point  is  dis¬ 
tributed  over  particles  of  different  sizes  and  fall  speeds,  V  at 
that  point  must  be  an  average  value.  Then  a  development  practically 
identical  to  that  in  Haurwitz  ^4,  p.  1 2qJ ,  for  example,  gives  the 
following  fundamental  continuity  equation  for  precipitation  content  M, 

(1) 

The  horizontal  winds  u  and  v  are  in  the  x-  and  y-  directions,  respec¬ 
tively,  and  the  vertical  air  speed  is  w  .  this  equation  governs  the 
distribution  of  M  already  formed  but  makes  no  allowance  for  the 
creation  of  M  (creation  of  M  is  considered  in  section  1,1.3  below). 

Assume  that  horizontal  and  local  time  changes  of  air  density  f 
are  relatively  small.  Then  the  air  motions  are  in  close  accord  with 
the  equation 
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Substitution  of  Eq.  (2)  into  Eq.  (l)  yields 
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The  derivatives  of  the  wind  appearing  in  Eq.  (l)  are  absent  from  Eq. 
(3);  the  density  tern  which  replaces  then  accounts  for  the  compressi¬ 
bility  of  the  atmosphere. 

1.1.2  A  continuity  equation  for  cloud  and  vapor 

A  tendency  to  supers&turation  in  the  atmosphere  is  rapidly  off¬ 
set  by  condensation  of  tiny  cloud  particles;  cloud  particles  in  clean 
air  evaporate  in  a  matter  of  seconds  when  the  relative  humidity  ie 
more  than  a  few  percent  below  saturation  £5,  7J.  these  properties  of 
cloud  suggest  our  definition  of  m:  m  is  the  cloud  density  minus 
the  saturation  vapor  density  plus  the  actual  vapor  density,  i.e. , 


m  ■*  cloud  density  +  Q  -  «s  -  2?-  q,  ,  U) 


where  Q  is  the  actual  vapor  density,  Q  is  the  saturation  vapor 

s 

density,  and  is  the  total  content  of  water  substance  excepting 

precipitation. 

A  continuity  equation  for  m  is  derived  by  first  considering 
that  for  ,  which  is  similar  to  Eq.  (3),  except  for  the  absence 

of  terms  like  M  -4%  and  V  : 
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Substitute  for  ^2  in  (5)  from  (4)  to  obtain 
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Assume  that  the  saturation  vapor  density  of  the  air  is  locally  steady 
and  horizontally  uniform,  i.e.,  that 

Even  though  this  is  rarely  exactly  true,  it  is  usual  for  the  local 
time  changes  and  horizontal  variations  of  Qg  to  be  quite  small  com¬ 
pared  with  the  vertical  variations  £cf.  assumptions  implicit  in  Eq. 
(2)J.  Then  Eq.  (6)  becomes 

it  ■  — 


Equation  (?)  is  significantly  unlike  Eq.  (3).  Since  m  has  no 
fall  speed  relative  to  the  air,  Eq.  (7)  does  not  contain  terms  in 
the  fall  speed  of  m;  also,  since  the  saturation  deficit  or  cloud 
content  changes  when  the  air  is  displaced  vertically,  Eq.  (7)  con¬ 
tains  terras  in  the  saturation  density  Qb  which  account  for  these 
changes.  The  sum  of  the  terras  in  Qs  in  (7)  we  call  the  generating 
function  G;  a  fair  approximation  to  G  in  the  real  troposphere  is 
a  linear  function  of  height. 

1.1.3  A  system  of  continuity  equations  for  vapor,  cloud,  and  preci¬ 
pitation  Including  cJ oud-precipitation  interactions 

The  modification  of  Eqs  (3)  and  (7)  to  model  the  bulk  effects 
of  interactions  among  the  three-dimensional  cloud  and  precipitation 
particles  as  well  as  their  distribution  by  wind,  is  based  on  elemen¬ 
tary  considerations  of  precipitation  physics.  The  first  product  of 
rising  motion  is  cloud.  Coalescence  of  cloud  particles  to  form 
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precipitation  is  apparently  favored  by  the  presence  of  a  broad  spec- 
trum  of  cloud  particle  sizes,  e.g.,  by  the  presence  of  large  hygro¬ 
scopic  salt  nuclei  in  the  subsaturated  atmosphere  .  Growth  of 
tiny  ice  particles  to  precipitation  size  may  also  occur  by  deposi¬ 
tion  when  liquid  and  ice  phases  coexist  at  subfreezing  temperatures 
^6,  7J .  Hither  of  these  processes  may  be  represented  by  a  term  of 
appropriate  magnitude  called  "autoconversion  of  cloud." 

Once  precipitation  particles  are  formed,  their  fall  speed  carries 
them  through  the  population  of  cloud  particles  and  facilitates  the 
rapid  growth  of  precipitation  by  collection  of  cloud  ^sj.  This  pro¬ 
cess  is  represented  by  a  term  whose  magnitude  must  increase  with 
the  content  of  cloud  and  precipitation;  this  term  called  "collection 
of  cloud",  contributes  equally  to  development  of  precipitation  and 
depletion  of  cloud. 

The  third  processes  of  particular  interest  is  evaporation.  The 
terms  in  Qg  in  Eq.  (7)  require  that  evaporation  of  cloud  in  down- 
drafts  is  just  rapid  enough  to  maintain  saturation.  The  evaporation 
of  precipitation  in  unsaturated  (cloud  free)  air  requires  an  additional 
term. 

Equations  (3)  and  (7)  are  now  rewritten  with  the  addition  of  terms 
representing  the  transfer  of  precipitation,  cloud,  and  vapor  from  one 
phase  to  another. 
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+  autocon.  of  cloud  +  collection  of  cld.  -  evap.  of  precip. 
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M  -  -«*■-  *5?  -  +  -  (9, 

-  autocon.  of  cloud  -  collection  of  cld.  +  evap.  of  precip. 

3h  Eq.  (9),  G  ■  »  where  Qg  is  the  satur¬ 

ation  vapor  density,  and  is  the  saturation  mixing  ratio.  Die 
autoconversion  and  collection  terms  in  (8)  and  (9)  are  non-zero  only 
when  m  ">  0,  i.e.,  when  there  is  cloud.  As  noted  above,  the  evapora¬ 
tion  term  is  non-zero  only  when  a  <0,  i.e.,  when  the  air  is  not 
saturated  with  vapor. 

Equations  (8)  and  (9)  can  be  reformulated  to  refer  to  mixing 
ratio  units  by  substituting  ^  m'  wherever  m  appears,  ^  M' 
wherever  M  appears,  and  0Q'  wherever  Q  appears.  After  simpli¬ 
fication,  there  is  obtained  * 


-■3F  -  "if  -(-  +  v)6f 

+  autocon.  of  cld.  +  cld.  collec.  -  evap.  of  precip. 


.u^al  _  v4al  .  w £sL  _ 
u  <?*  w  a* 

-  autocon.  of  cld.  -  cld.  collec.  +  evap.  of  precip. 


(9a) 


where  the  primed  quantities  are  in  mixing  ratio  units. 

The  absence  of  the  density  term  in  (9a)  is  due  to  the  conserva¬ 
tive  nature,  following  the  motion  of  air  particles,  of  the  mixing  ratio 
of  water  vapor  plus  cloud.  Although  (9a)  is  therefore  slightly  simpler 
than  Eq.  (9),  our  studies  use  Eqs.  (8)  and  (9)  in  density  units  be¬ 
cause  radar  reflectivity  characteristics,  visual  appearance  of  clouds, 
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and  certain  physical  effects  are  more  easily  understood  in  these  units 


The  published  works  ^2,  3,  4j  which  form  a  principal  foundation 
to  the  studies  of  this  contract,  treat  a  cloud-free  model  atmosphere 
in  which  condensate  appears  instantly  as  precipitation  of  fall  speed 


V.  Equations  (8)  and  (9)  can  be  reduced  to  the  basic  equation  of  the 
earlier  model.  Thus,  the  addition  of  Eqs.  (8)  and  (9)  yields 


st 


_  u  +  M)  _  v  „  Msl±M  _  -2.  „v 

u  ax  v  St  st 

+  (m  +  M)  w  +  wG  . 


If  M  +  ■  is  denoted  by  a  new  variable  M*  and  if  it  is  assumed  that 
m  *=  0  when  M*  >  0  and  that  M  and  V  are  zero  when  0,  Eq. 

( 1 )  becomes 

am  .  . .  .  T  4a; . .  ^*e .  +  mv -aaf  +  -) 

at  dx  ST  Sz  Sz  Szs 

(11) 

(m  “  0  when  M*  >  0 
)M  -  V  -  0  when  M*  C  0 


This  is  the  three-dimensional  formulation  of  the  model  system  that  has 
been  treated  extensively  in  the  earlier  work.  M*  when  negative  is 
the  saturation  deficit  which  follows  the  motion  of  the  air;  when  posi¬ 
tive,  M*  is  the  precipitation  content  which  falls  at  speed  V  rela¬ 
tive  to  the  air.  Equation  (ll)  models  a  system  wherein  cloud  does  not 
exist  because  condensate  appears  instantly  as  precipitation. 
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1 .3  Steady-state  Vertical  Distribution  of  Cloud  in  a  Precipitation- 

free  Compressible  Atmosphere 

Section  1.3  in  Report  No.  2  discusses  the  steady-state  vertical 
profile  of  cloud  in  a  precipitation-free  compressible  atmosphere  as 
modeled  by  the  equation 

-If  -  .^ff  +  o.  <-*> 

The  compensation  level  is  defined  as  that  height  where  m ft  *  -  0, 
i.e.,  where  the  tendency  of  the  generating  function  to  increase  the 
density  of  condensate  is  exactly  compensated  by  the  decreasing  density 
of  ascending  air.  Above  the  compensation  level,  cloud  density  decreases 
with  ascent  and  increases  with  descent.  In  Report  No.  2,  a  function 
linear  with  height  is  assumed  for  G  ,  and  the  corresponding  solutions 
of  Eq.  (12)  are  presented  in  Fig.  4  of  that  report.  After  Report  No.  2 
was  written,  a  more  accurate  profile  of  the  steady-state  cloud  content 
was  derived.  Since  no  relatively  simple  formula  has  been  found  to  approx 
imate  the  natural  generating  function  at  high  levels  (low  temperatures) 
with  the  desired  accuracy,  the  revised  cloud  profiles  have  been  computed 
with  the  aid  of  the  Smithsonian  Meteorological  Tables  jV,  Tables  71,  78, 
and  108  . 

The  results  for  an  atmosphere  of  constant  wet-bulb  potential  tem¬ 
perature  295°  are  shown  in  Fig.  1.  Note  that  in  the  upper  troposphere, 
(ice)  cloud  amounts  less  than  0.1  gm  may  be  associated  with  anomalous 
behaviour,  i.e.,  decreasing  cloud  density  with  ascending  motion,  and  con¬ 
versely,  Accurate  measurements  of  water  content  and  temperature  in  dense 
tropical  cirrus  would  be  of  special  in.arest  in  connection  with  the  appli 
cation  of  this  theory  to  problems  concerning  the  development  and  persis¬ 
tence  of  high  cloud. 


Fig.  1  Steady-state  cloud  profiles  in  a  non- precipitating  at¬ 
mosphere  with  wet  bulb  potential  temperature  -  22°C. 
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Height,  km  Height, 


The  data  represented  by  Fig.  1  vary  only  slowly  with  height  for 
constant  temperature;  this  figure  can  therefore  be  applied  to  various 
elevations  where  the  same  temperatures  occur  in  temperate  and  arctic 
regions . 

1 . L  A  Model  of  the  Chset  of  Precipitation 

lhe  cloud  conversion  and  accretion  terras  in  (8)  and  (9)  have 
been  discussed  in  the  previous  reports  on  this  contract.  The  natural 
processes  are  modeled  as  follows: 


Cloud  autoconversion  -  k1 (m  -  a) 
Accretion  ■  kg&iQ*12^ 


(k,  •*  0  when  m  <  a)  .  . 

/  v  (13) 

(k2  “  0  when  m  <  0) 


In  (13),  k1  and  a  are  parameters  whose  magnitudes  are  selected  to 
model  spontaneous  conversion  of  cloud  at  rates  in  accordance  with 
theories  and  observations,  kj  •  7  I  10”\  nQ  defines  the  number  of 
small  drops  in  an  assumed  Marshall-Palmer  distribution  of  precipitation 
particles,  and  E  is  the  efficiency  with  which  precipitation  particles 
collect  cloud.  Units  are  gns ,  meters,  and  seconds  throughout. 

In  section  1.5  of  Report  No.  2  the  onset  time  of  precipitation  in 
the  incompressible  one-dimensional  model  version  of  Eqs.  (8)  and  (9) 
is  related  to  the  magnitude  of  k^ .  For  values  of  the  other  parameters 
which  are  considered  to  be  reasonable  approximations  to  natural  values, 
the  onset  time  of  precipitation  at  the  ground  was  found  to  increase 
only  about  10  minutes  with  a  change  of  k^  from  10”1  sec-1  to  10”'* 
sec”^ .  In  order  to  understand  better  how  the  cloud  conversion  and 
accretion  coefficients  influence  the  onset  of  precipitation,  exact 
solutions  of  a  simple  model  which  embraces  only  the  initiation  problem 
are  studied. 
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Consider  the  terns  (13 )  as  they  operate  in  a  deep  Model  cloud 
of  unifora  water  content  without  vertical  air  motion.  In  this  model, 
cloud  changes  to  precipitation  everywhere  at  the  same  rate,  the  ad- 
vection  terms  vanish,  and  the  equations  describing  the  amounts  of 
cloud  and  precipitation  are 

«  -  -dt  "  >=,(--•)  +  W'25  -r875  (14) 

m  0  when  a  <a; 
k1  >  0  when  a  >  a; 
k2  >  0  always  . 

Since  there  is  no  vertical  motion  and  therefore  no  condensation,  the 
condition 

M  +  m  “  m(t  ■  0)  ■  mQ  ,  0  ^  t  £:  (l 5 ) 

where  Bq  is  the  initial  cloud  amount,  must  also  be  satisfied.  If 
the  exponent  of  M  in  ( 14 )  is  changed  from  *f  0.875  to  +  1.0  ,  exact 
integration  is  grately  facilitated  without  changing  essential  proper¬ 
ties  of  the  formulation.  Then  substitution  of  (l 5 )  into  (14)  with 
altered  exponent  yields 

dt  "  '  dt  *  k1  (m  -  a)  +  C^Osq-m),  (16) 

where  C2  -  k2  EhQ*125. 

The  interval  t^  required  for  the  cloud  amount  to  decrease  from 
the  amount  nig  to  the  amount  a  ,  when  the  autoconversion  term  ceases 
to  be  a  factor,  is  a  convenient  measure  of  the  time  of  onset  of  precipi¬ 
tation  : 
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■•-f 

Jnr 


C2"2  ‘  (k1  +  C2n0)  B  +  V 


(17) 


The  Integral  of  the  right  hand  side  depends  on  the  size  of  the  dis¬ 
criminate  j0  of  the  denominator 


0 


k,2  +  “i0*  *  iV  -  “iV 


(18) 


That  0  >  0  can  be  shown  in  the  following  way.  Note  first  that 
k, ,  C^,  mQ,  and  a  are  all  positive.  Therefore,  the  first  three  terns 
on  the  right  hand  side  are  positive.  Set  mn  *  a;  if  y0>O  when 


a,  the  same  is  certainly  true  when  m0*>a  .  Thus 
£)  >  kj2  +  2k,C2a  +  C2a2  -  Ak^a 


and 


0 


>  (k,  -  c2.r 


(19) 


(20) 


(fcvioualy  the  minimum  value  of  the  expression  on  the  right  hand  side 
is  zero  and  the  discriminant  is  therefore  positive.  Then 

l-i 

|i 
li 


t  » 
a 


(ki  *  cz"o)2  -  “li*]' 

ta  2C2»  -  <kl  C2»0>  ~  [  <k1  *  i»0 )2  -  “l  C2»] ' 
20^-  (k,  +  C^q)  +  [(k,  +  C^q)2  -  4k, C2a] f 


.(21) 


When  8q  *=  2a,  this  can  be  written 


ta  -  J  In 


(22) 


where 
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(k,  +  4C22a2)*  , 


ip  2*2 
4C0  a  i 

-1  -  (1  +  — ^-r-P 

k 

_ K1 

,  c  2,2 
4u_  *  1 

-1  +  (1  +  -S-s-)* 
k1 


(3) 


Note  that  S  and  j*/  for  specified  a  depend  only  on  the  ratio 
C2a/k1  ;  all  the  terms  are  dependent  on  k^  alone  when  the  term  C,,a 
is  fixed. 


The  variation  of  t  with  C  a  for  m~  **  2a  and  for  various 
a  <c  u 

values  of  the  cloud  conversion  coefficient  k1  is  shown  in  Fig.  2. 
Note  that  when  k^  >>  C^a  ,  the  time  tft  is  a  very  slowly  varying 
function  of  C^a.  And  when  k1  <<  C^a  ,  t&  is  a  verjt  slowly  varying 
function  of  k^  .  In  other  words,  if  the  coalescence  and  aggregation 
processes  among  the  many  tiny  cloud  drops  are  relatively  rapid,  the 
efficiency  with  which  the  larger  raindrops  collect  cloud  is  of  little 
consequence  to  the  rate  of  depletion  of  cloud:  and  if  the  collection 
of  cloud  by  the  raindrops  is  relatively  rapid,  then  the  magnitude  of 
the  rate  of  cloud  conversion  is  not  critical,  so  long  as  it  is  more 
than  zero. 
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1 
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0.125 


Fig.  2  Variation  of  the  onset  parameter  ta  with 
the  conversion  parameter  k-|  and  with  accretion. 
It  is  assumed  that  idq  **  2a. 


Fig.  3  Detailed  development  with  time  of  precipi¬ 
tation  in  the  model  defined  by  Eq.  (16)  with 
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_-3  _2 

In  Fig.  2,  the  region  10  J 4.  C^a  <  10  holds  the  greatest 
interest  because  it  is  here  that  most  empirical  data  are  represented. 
For  example,  the  heavy  horizontal  line  indicates  the  following  parame¬ 
ter  combination: 

—3  -3  7 

a  =  1  gm  m  ,  Hq  =  2  gm  m  ,  E  *=  1,  and  Hq  *=  10  . 

—3  —2  7  -A 

If  a  ■=  0.5  gm  a  ,  Bq  *  1  gm  a  ,E  =  .5,  and  nQ  =  10  m  4  then 

C^a  =  1.25  X  10"3.  The  values  of  in  the  range  10“^  kj  6  10-1 

—1  -3  -2 

sec  for  the  range  10  C^a  £  10  are  associated  with  onset 

parameters  t  ranging  from  about  one  minute  to  about  on  hour. 

When  mis  a,  k1  =0  and  the  solution  of  (16)  is  simply 


This  equation  has  contributed  to  the  derivation  of  the  curves  in  Fig. 

3  which  illustrate  in  detail  the  time  dependence  of  the  cloud  content 
for  three  values  of  C^a/k^  .  Curves  of  the  precipitation  content  M 
mirror  those  shown;  they  start  at  M  =  0  and  asymptotically  approach 
the  value  Hq  as  the  time  increases. 

1.5  Relationships  Between  Microphysical  Parameters  and  Transient 
Features  of  Solutions  of  the  Complete  Model 

The  one-dimensional  forms  of  Eqs.  (8)  and  (9)  for  an  incompressi¬ 
ble  atmosphere  have  been  solved  by  an  IBM  7094  programmed  to  use  the 
finite-difference  equations  discussed  in  Report  No.  1.  In  each  case 
discussed  here,  the  initial  condition  is  assumed  to  be  a  saturated 
atmosphere  without  cloud  or  precipitation.  The  updraft 
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(a  -  z  /H),  where  w  m  0.5  n/aec  and  H  ,  the  top  of  the  updraft 
n&x 

column,  is  6  km  in  all  cases.  The  generating  function  G  -  3  x  10  - 

_7 

3  X  10  a  ,  is  a  fair  approximation  to  the  lower  troposphere  in  the 
tropics.  This  section  discusses  the  role  of  the  model  microphysical 
parameters  in  shaping  transient  features  of  the  solutions. 

1-5.1  The  onset  parameter  a 

The  model  provides  that  conversion  of  cloud  to  precipitation  does 
not  start  until  the  cloud  content  produced  in  saturated  updrafts  has 
exceeded  the  magnitude  a  .  Thereafter,  cloud  responds  to  its  accre¬ 
tion  by  relatively  large  precipitation  particles,  in  addition  to  ad- 
vection  and  condensation  processes.  Often,  the  amount  of  cloud  which 
can  coexist  in  equilibrium  with  precipitation  is  less  than  the  amount 
required  to  initiate  precipitation  in  the  model.  When  this  is  the 
case,  the  time  dependent  solutions  for  cloud  and  precipitation  are 
characterized  by  an  initial  pulse  whose  height  is  about  proportional 
to  the  difference  between  the  equilibrium  precipitation  content  and 
the  magnitude  of  a  . 

A  mathematical  approach  to  this  problem  can  be  illustrated  for 
a  place  on  the  cloud  profile  where  the  cloud  content  is  at  or  near  its 
maximum  and  the  advection  of  cloud  is  therefore  small,  and  where  the 
steady-state  precipitation  rate  is  nearly  eoual  to  the  condensation 
rate  integrated  over  all  higher  elevations.*  Then 

*  i.e.,  the  equilibrium  cloud  content  is  small  compared  to  its  steady- 
state  value  in  the  absence  of  precipitation,  and  updrafts  are  small 
comnared  to'  precipitation  fall  speeds  (cf.  Report  No.  2,  Section  1.2). 
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K 


-  M(V  +  w)  » 


wGdz  , 


(25) 


where  R  is  the  precipitation  rate  at  height  z  j  V  ,  the  fall  velo- 
city  of  M  ,  discussed  in  the  earlier  reports,  is  V= -38.6  nQ  ’  M’  , 
and  H  is  the  top  of  the  saturated  updraft  column.  An  equilibrium 
between  condensation  and  accretion  processes  is  represented  by  the  equa¬ 
tion 

0  =  wG  -  k„  En  *125  mM  .  (26) 


By  carrying  out  the  integration  in  (25)  for  a  parabolic  w-  profile  and 
the  generating  function  G  =  A  +  Bz  ,  and  making  other  obvious  sub¬ 
stitutions,  Eqs.  (25)  and  (26)  can  be  combined  to  obtain  the  following 
expression  of  the  maximum  cloud  amount  which  can  coexist  in  equi¬ 

librium  with  precipitation: 


in 


3.35  X  104  /*maxx 
E 


2/9 


"ep  E 

In  Eq.  (27),  i  and  Y 

4*  =  AHz  -  (A  -  BH)  z2  -  Bz3 


, f  -7/9 

4 

\i/  ' 

(27) 


are, 


(28) 


*VK  “  fAH(H2-  z2)/^-[(A-  HBHH3-  z3)/^-[b(H^-  z4)/a]  . 


It  is  interesting  that  greatly  simplifying  assumptions  still  lead 
to  quite  complicated  expressions.  The  magnitude  of  the  initial  pulse 
of  cloud,  which  is  manifested  in  the  time  variations  of  precipitation 
at  the  surface,  can  be  estimated  by  comparing  the  magnitude  of  m 
given  by  (27)  with  the  magnitude  of  a  .  Plots  of  Eq.  (27)  will  be 
presented  in  the  next  report. 
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1.1.2  The  cloud  conversion  parameter  k. 


In  the  model,  the  parameter  is  the  rate  of  autocomrersion 
to  precipitation  of  cloud  content  in  excess  of  the  magnitude  a  . 

When  kj  is  large,  the  approach  to  a  steady-state  condition  is  ex¬ 
pedited  after  m  exceeds  a  ;  when  is  small,  the  cloud  content 
continues  to  increase  in  the  updrafts  after  precipitation  starts 
(especially  if  C^a  is  also  small).  Thus  the  height  of  a  pulse¬ 
shaped  transient  near  the  start  of  precipitation  is  enhanced  as  the 
value  of  kj  decreases.  This  is  illustrated  for  three  values  of 
kj  in  Fig.  4a.  Nearly  the  same  steady  state  is  ultimately  approached 
in  all  three  cases  because  the  cloud  content  in  equilibrium  with  pre¬ 
cipitation  for  the  indicated  choice  of  parameters  is  less  than  the 
amount  a  at  most  levels,  and  the  magnitude  of  is  therefore 
rather  unimportant  after  the  steady-state  condition  in  these  cases 
has  been  established.  (See  Fig.  5b,  section  1.6.2.) 

1.5.3  The  collection  efficiency  E 

The  model  permits  study  of  the  role  in  the  precipitation  process 
of  the  efficiency  with  which  precipitation  collects  cloud.  Figure  4b 

I 

shows  the  development  of  precipitation  at  the  ground  in  a  case  where 
the  collection  efficiency  E  is  unity  and  in  a  case  where  E  =  0. 

The  steady-state  precipitation  rate  is  lower  in  the  latter  case,  be¬ 
cause  more  condensate  remains  as  cloud  aloft  and  is  spread  horizon¬ 
tally  by  high  level  divergence  of  the  wind.  When  E  is  greater  than 
zero,  two  processes  can  contribute  to  precipitation  formation;  when 
E  is  relatively  very  small  or  zero,  the  autoconversion  process  is 
the  only  important  one. 
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Fig.  4a  Development  of  precipitation  at 
the  ground  for  various  magnitudes  of  the 
conversion  parameter  . 


Fig.  4b  Development  of  precipitation  at 
the  ground  for  extrema  of  the  collection 
efficiency  E . 


Fig.  4c  Development  of  precipitation  at 
the  ground  for  three  values  of  n_  .  The 
value  n_*  ictfm-i  is  typical  of  natural 
rains. 
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1.5.4  The  precipitation  drop  tin  Perimeter 

For  at  least  two  reasons,  the  parameter  tiq  affects  the  cloud 
and  precipitation  distributions.  First,  the  rate  of  cloud  collection 
for  fixed  M  and  E  increases  slowly  as  nQ  increases.  Second,  as 
n^  increases,  the  number  of  small  drops  increases,  and  therefore  the 
average  terminal  velocity  of  drops  decreases.  However,  both  the  collec¬ 
tion  and  fall  speed  equations  involve  the  eighth  root  of  nQ  ,  and 
order  of  magnitude  changes  of  this  distribution  shape  parameter  are 
necessary  to  alter  cloud  and  precipitation  distributions  significantly. 

Figure  4c  shows  the  time  variation  of  precipitation  at  the  ground 
accompanying  three  values  of  n^.  The  middle  curve  is  associated  with 
an  nQ  typical  of  natural  rain.  When  nQ  is  very  small,  particles 

fall  rapidly  and  precipitation  starts  earlier  at  the  ground.  When  nQ 
10  -4 

is  as  large  as  10  m  ,  however,  the  average  particles  are  only  the 
size  of  drizzle;  their  very  slow  falling  speed  is  associated  with  both 
the  delayed  arrival  of  steady  conditions,  and  the  relatively  great 
rate  of  steady  precipitation  of  this  case.  The  latter  is  associated 
with  the  air's  horizontal  convergence,  implicit  in  the  equations,  which 

I 

produces  a  progressively  greater  imbalance  between  condensation  and 
precipitation  rates  at  updraft  centers  as  the  ratio  of  fall  speeds  to 
updrafts  approaches  unity.  As  discussed  previously  (section  1 .2  in 
Heport  No.  2)  an  excess  or  deficit  of  precipitation  at  updraft  centers 
as  compared  to  condensation  there,  is  compensated  at  locations  away 

*'•  A* 

from  the  center. 
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1 .6  Relationships  Between  Hicrophvsical  Parameters  and  Steady-state 

Profilea  of  Cloud  and  Precipitation 

Hula  section  discusses  the  role  of  the  microphyslcal  parameters 
in  shaping  steady-state  solutions  of  the  differential  equations.  The 
initial  conditions  and  the  parameters  ,  H  and  0  are  the 

same  as  noted  in  section  1.5.  In  section  1.5  the  discussion  of  tran¬ 
sients  is  illustrated  by  graphs  of  the  development  with  time  of  the 
precipitation  rate  at  the  ground.  In  this  section,  the  figures  show 
steady-state  vertical  profiles  of  cloud  and  precipitation. 

1.6.1  The  onset  parameter  a 

When  the  onset  parameter  a  is  quite  large,  the  cloud  amount 
is  determined  solely  by  a  balance  between  condensation,  vertical 
advection,  and  accretion  processes.  When  a  is  small,  however,  the 
autoconversion  and  accretion  processes  both  contribute  to  depletion 
of  cloud,  and  the  steady-state  cloud  content  is  therefore  reduced. 

The  effect  of  varying  a  between  particular  values  is  shown 
in  Fig.  5a.  When  a  is  large,  the  cloud  content  in  the  upper  atmos¬ 
phere  tends  to  be  large  because  it  is  here  that  precipitation  ie 
light  and  relatively  ineffective  in  removing  cloud.  some  of  the 
large  amounts  of  cloud  aloft  are  lost  to  precipitation  at  the  up¬ 
draft  center  by  virtue  of  horizontal  divergence  of  the  wind  £cf, 

Eq.  (7)  in  Report  No.  2J,  and  the  precipitation  at  the  ground  is 
therefore  slightly  reduced  in  this  case. 

Note  that  in  both  cases  the  cloud  amount  below  3.5  km  is  less 
than  a  and  therefore  locally  regulated  by  the  accretion,  conden¬ 
sation,  and  advection. 
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1.6.2  The  cloud  conversion  parameter  kj 


Comparison  of  the  steady-state  profiles  in  Figs.  5a  and  5b 
shows  that  the  roles  of  k^  and  a  are  similar.  Indeed,  the  cloud 
and  precipitation  profile  in  Figs.  5b  for  m  10  a  *=  0.5  is 
virtually  identical  to  the  profile  for  kj  ”  10  a  =  2.0  in  Fig. 
5a.  The  similarity  is  not  surprising,  since  autoconversion  is  made 
small  either  by  increasing  a  or  decreasing  k1  . 

1.6.3  The  collection,  efficiency  E 

The  role  of  the  collection  efficiency  E  is  illustrated  by  Fig. 
5c.*  Whan  the  collection  efficiency  is  zero,  the  accretion  process 
no  longer  contributes  to  depletion  of  cloud,  and  the  cloud  content 
rises.  The  curve  shown  for  E  =  0  applies  to  a  balance  between  con¬ 
densation,  advection,  and  cloud  autoconversion  processes. 

Note  that  the  change  of  the  cloud  profile  for  the  change  of  E 
from  1  to  0  is  much  greater  than  the  accompanying  change  of  the  pre¬ 
cipitation  profile;  this  is  true  for  the  particular  values  of  k^ 
and  y  used  because  k^  alone  is  large  enough  in  this  case  to 
permit  most  of  the  cloud  to  convert  to  precipitation  before  being 
spread  by  horizontal  divergence  at  high  levels.  In  other  words,  the 
conversion  time  constant  1/k  is  much  smaller  (l/k  =  1000  secs.) 
than  the  time  that  it  takes  an  air  parcel  to  ascend  from  low  to  high 

levels  in  the  column  (H/w  =  12,000  secs.). 

max 

By  setting  =  0  in  the  one-dimensional  incompressi¬ 

ble  form  of  Eq.  (9i  one  finds  that  when  in  the  absence  of  accretion, 
the  cloud  profile  has  a  maximum,  that  maximum  satisfies  the  equation 
*  The  curve  for  E  «=  1'  in  this  figure  also  appears  in  Figs.  5a  and  5b, 
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Fig.  5a  Steady-state  profiles  of  cloud 
and  precipitation  content  for  two  values 
of  the  onset  parameter  a. 


Fig.  5b  Steady-state  cloud  and  precipi¬ 
tation  content  profiles  for  three  values 
of  the  cloud  conversion  parameter  . 


Fig.  5c  Steady-state  cloud  and  precipi¬ 
tation  content  profiles  for  extrema  of 
the  collection  efficiency  E . 


Fig,  5d  Steady-state  cloud  and  precipi¬ 
tation  content  profiles  for  three  magni¬ 
tudes  of  T\q  . 
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(29) 


ataax 


If  the  location  of  were  fixed,  its  magnitude  would  rise  in¬ 

definitely  as  decreases.  Equation  (29)  is  indeterminate,  how¬ 
ever,  unless  the  location  of  m  and  the  magnitude  there  of  wG 

max 

are  ascertained.  The  maximum  numerical  value  of  m  at  any  particu¬ 
lar  point  can  not  be  larger  than  the  steady-state  cloud  content  at 
that  point  in  a  nonprecipitating  atmosphere.  As  becomes  smaller, 
the  m-  profile  in  the  absence  of  accretion  more  nearly  resembles  the 
steady-state  profile  defined  by  the  equation 


this  equation  has  m  increasing  with  z  to  the  top  of  the  updraft 
column. 

When  1  /k1  >  ">  ,  the  cloud  profile  is  accurately  defined 

by  Eq.  (30)*  and  condensation  is  spread  horizontally  by  high  level 
divergence  instead  of  producing  rain  at  places  beneath  the  updraft 
column. 

1.6.4  The  precipitation  drop  size  parameter  ng 

Much  of  the  explanation  of  the  curves  shown  in  Fig.  5d  is  given 
in  section  1.5.4.  When  nQ  is  large,  small  slow  falling  drops  are 
numerous  and  accretion  of  cloud  is  enhanced.  Therefore,  cloud  amounts 


*  In  the  compressible  case,  the  steady-state  m-  profiles  in  the  ab¬ 
sence  of  accretion  and  as  kj  approaches  0  ,  approach  the  curves 
discussed  in  section  1.3  of  this  report. 
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in  this  csss  are  small.  The  precipitation  content  M  is  increased 
when  nQ  is  large  and  fall  speeds  small  because  there  is  then  more 
time  for  the  precipitation  to  develop  from  condensation  in  the  up¬ 
drafts.  And  the  enhancement  of  the  precipitation  rate  MV  at  up¬ 
draft  centers  when  fall  speeds  are  small  can  be  explained  by  the 
divergence  term  in  Eq.  (?)»  Report  No.  1,  whose  contribution  to  the 
surface  precipitation  rate  increases  m  the  ratio  of  precipitation 
fall  speeds  to  maximum  updrafts  approaches  unity. 

1.7  Computer  Program  for  Two-dimensional  Rectilinear  and  Radially 
Synetric  Models 

Dr.  Newburg  is  preparing  specifications  for  a  computer  program 
to  define  cloud  and  precipitation  distributions  and  budget  parameters 
In  two-dimensional  rectilinear  and  radially  symmetric  model  wind  fields. 
An  explicit ,  conditionally  stable  finite  difference  approximation 
to  the  partial  differential  equations  will  be  used  in  this  program. 
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2.0  CONCLUSIONS 

At  heights  greater  than  12  km  in  the  tropical  atmosphere,  cloud 
amounts  smaller  than  .05  gm/rn^  are  affected  only  slightly  by  verti¬ 
cal  displacements  of  the  order  of  a  kilometer.  Ibis  fact  may  explain 
relatively  great  persistence  of  cirrus  cloudiness. 

The  cloud  and  precipitation  model  illustrates  that  when  more 
than  one  process  influences  precipitation  production,  the  onset  of 
precipitation  is  controlled  only  within  certain  limits  by  alteration 
of  one  of  the  processes.  Where  two  processes  for  precipitation  forma¬ 
tion  exist,  simultaneous  changes  in  the  same  direction  of  the  magni¬ 
tude  of  both  are  required  if  the  changes  are  to  be  linearly  reflected 
in  the  time  of  onset  and  amount  of  precipitation. 

The  discussion  of  sections  1.5  and  1.6  show  that  the  complete 
kinamatical  equations  with  specified  updrafts  is  a  model  system  whose 
reaction  to  changes  of  the  microohysical  processes  is  generally  in  a 
direction  to  minimize  the  effects  of  the  changes  on  the  surface  preci¬ 
pitation  rate,  Thus,  an  increase  of  cloud  collection  efficiency 
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is  associated  more  with  decreased  cloud  amounts  than  with  increased 
precipitation.  Although  the  fall  speed  of  precipitation  particles 
is  closely  associated  with  the  distribution  of  precipitation  in  the 
neighborhood  of  an  updraft  center,  the  total  amount  of  condensation 
and  precipitation  is  most  closely  related  to  the  strength,  areal  ex¬ 
tent,  and  duration  of  updrafts. 

The  work  should  be  of  value  for  interpreting  measurements  of 
the  development  of  cloud  and  precipitation  distributions  in  terms  of 
the  accompanying  wind  fields  and  microphysical  processes.  It  may 
also  stimulate  weather  modification  experiments  designed  to  alter 
precipitation  and  cloud  distributions. 

The  combination  of  thermohydrodynamic  considerations  with  the 
kinematic  theory  would  be  an  important  achievement. 

3.0  PROGRAM  FOR  THE  NEXT  INTERVAL 

The  following  work  and  study  has  been  planned  for  the  next  in¬ 
terval: 

1 )  Study  of  the  development  of  a  model  cloud-precipitation  sys¬ 
tem  from  unsaturated  initial  conditions. 

2)  Examination  of  equations  which  relate  the  rate  of  evaporation 
of  precipitation  beneath  a  convective  cloud  to  the  saturation  deficit, 
the  rate  of  precipitation,  and  the  vertical  velocity  profile. 

3)  Preparation  of  a  computer  program  to  derive  cloud  and  preci¬ 
pitation  distributions  and  budget  parameters  in  two-dimensional  rec¬ 
tilinear  and  radially  symmetric  model  wind  fields. 

4)  Preparation  of  the  one-dimensional  results  for  submission  to 
the  Journal  of  Atmospheric  Sciences. 


26 


L.O  PARTICIPATICM  OF  PERSONNEL 
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